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1.0  EXECUTIVE  SUMMARY 


A  preliminary  design  concept  has  been  developed  for  a  holographic  instrument  to 
measure  scattering  and  beam  transmission  optical  properties  of  ocean  waters.  The 
focus  of  this  design  process  is  an  instrument  which  can  be  developed  into  an 
expendable  unit  with  size  constraints  imposed  by  the  expendable  batho- 
thermograph/airbome  expendable  batho-thermograph  (XBT/AXBT)  technology.  The 
instrument  design  has  been  separated  into  four  measurement  functions:  (1)  beam 
transmission,  (2)  small  angle  scattering,  (3)  large  angle  scattering,  and  (4) 
backscatter.  The  present  research  has  focussed  on  developing  and  demonstrating  a 
holographic  design  to  measure  small  angle  scattering.  This  instrument  measures  the 
modulation  transfer  function  (MTF)  of  the  particles  suspended  in  the  water.  The 
MTF  can  be  inverted  with  a  Fourier-Bessel  transform  procedure  to  estimate  the 
volume  scattering  function  at  small  angles.  A  series  of  simulation  codes  were 
developed  to  investigate  instrument  performance.  In  a  laboratory  demonstration  a 
hologram  was  constructed  to  focus  a  sinusoidal  pattern  scattered  by  polystyrene 
spheres  onto  a  linear  detector  array.  The  laboratory  results  showed  loss  of  contrast 
with  increasing  spatial  frequency  for  particle  sizes  greater  than  1 .0  micron  diameter. 
The  laboratory  observations  were  consistent  with  the  model  simulations  both  for  the 
small  (0.49/im)  and  large  (25.0^m)  particles.  The  simulation  model  results  for  a 
25cm  instrument  demonstrated  the  capability  to  measure  contrast  loss  at  selected 
spatial  frequencies  and  used  that  information  to  derive  the  volume  scattering  function 
for  angles  from  0.005°  to  0.5°.  Simulation  of  the  laboratory  measurements  for  25 
micron  particles  were  consistent  with  unsealed  laboratory  measurements. 
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2.0  INTRODUCTION 


The  objectives  of  this  research  are  (1)  to  develop  a  concept  for  an  expendable 
optical  sensor  for  ocean  transmission  and  scattering  measurements  using  Holographic 
Optical  Elements  (HOEs),  (2)  to  evaluate  performance  by  simulating  the  instrument 
measurement  properties,  and  (3)  to  demonstrate  the  performance  of  the  HOE  in  the 
laboratory,  using  a  known  particle  size  distribution  and  scattering  functions  developed 
from  Mie  calculations.  The  first  section  of  the  report  describes  the  overall  instrument 
preliminary  design  which  has  been  partitioned  into  four  measurement  functions:  beam 
transmission,  small  angle  scattering  (0<1O"),  large  angle  scattering  (10°<^<170^, 
and  back  scattering  near  180°  (17O°<0<  180°).  For  the  present  study  the  research 
was  focussed  on  developing  a  holographic  small  angle  scattering  meter.  Following  a 
discussion  of  the  proposed  instrument  to  measure  ocean  optics  is  an  outline  of  the 
theory  for  small  angle  scattering  measurements.  The  next  section  is  a  description  of 
the  laboratory  setup.  Following  are  sections  describing  laboratory  and  model 
simulation  results.  The  last  section  of  the  report  provides  guidelines  and 
recommendations  for  development  of  a  integrated  holographic  instrument. 

2.1  OCEAN  OPTICAL  INSTRUMENT  OVERVIEW 

The  proposed  expendable  instrument  is  being  developed  for  ship  and  survey 
aircraft  deployment  (see  Figure  1).  The  airborne  version  will  be  launchable  from  a 
standard  sonobuoy  tube.  The  submersible  instrument  will  measure  inherent  optical 
properties  of  the  ocean  which  are  independent  of  incident  illumination  at  the  surface. 
The  principal  inherent  properties  are  absorption  (a),  beam  attenuation  (c),  and  the 
volume  scattering  function  a{6).  These  properties  together  specify  radiative  transfer 
in  the  oceanic  medium  and  are  interrelated  by  the  following  simple  expressions. 


OPERATIONAL  SEQUENCE 


1 )  The  Unit  is  Deployed  by  Hand  Over  the  Side  of  Any  Ship  Without  a 
Launcher  or  From  Aircraft  Deployed  From  a  Sonobuoy  Launch  Tube 

2)  The  Buoy  is  Activated  Electronically  Upon  Contact  With  the  Water 

3)  A  Carbon  Dioxide  Cartridge  Inflates  a  Flotation  Bag  of  the  Surface 
Module  and  Erect  and  Antenna  Within  the  Bag 

4)  The  Surface  Buoy  Drifts  Away  From  the  Ship  Before  Releasing  the 
Optical  Instrument  Thereby  Avoiding  Ship  Shadow  and  Wake 
Interference 

5)  The  Optical  Data  Obtained  by  the  Descending  Probe  is  Converted  to 
Frequency  Data  and  Telemetered  Via  Signal  Wire  to  the  Surface  Buoy 
Where  it  is  Transmitted  to  the  Ship  or  Aircraft.  The  Raw  Data  is 
Received  by  a  PC  Data  Logger  and  Displayed  in  Real  Time  as  the  Probe 
Descends 

6)  After  the  Probe  has  Completed  Its  Decent,  the  RF  Transmitter  is 
Disabled  and  the  Surface  Buoy  is  Scuttled  (Hannon,  1989) 


Figure  1.  Expendable  Instrument  Deployment 
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where  b  is  the  total  scattering  coefficient.  The  current  instrument  design  concept  has 
the  capability  to  measure  beam  attenuation  and  the  volume  scattering  function.  Use 
of  holographic  optics  can  meet  instrument  optical  performance  requirements  and 
provide  integration  of  complex  optical  components  into  a  compact  space  with  low 
production  cost. 

There  is  Navy  need  for  extensive  knowledge  of  the  optical  properties  of  the 
oceans  to  predict  performance  of  optical  systems  operating  in  the  ocean  environment 
and  to  expand  remote  observations  of  i  ocean  environmental  parameters.  In 
addition,  there  is  a  science  objective  to  obtain  optical  closure  between  attenuation, 
scattering,  and  absorption.  Oceanic  volume  scattering  functions  are  highly  anisotropic 
and  characteristically  have  a  large  forward  diffraction  peak.  Presently  there  are  very 
few  instruments  capable  of  small  angle  scattering  or  absorption  measurements.  Thus, 
development  of  a  small  angle  scattering  measurement  technique  and  instrument  is 
important  to  both  NAVY  science  and  operational  systems  objectives  for  ocean  optics. 
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2.2  HOLOGRAPHIC  OPTICAL  SYSTEMS 


The  need  for  complex  optical  systems  in  a  compact  expendable  ocean  instrument 
can  be  met  using  holographic  optical  elements.  A  hologram  is  a  coded  diffraction 
pattern  on  a  film  plate.  When  a  hologram  is  illuminated  with  coherent  light  of  the 
same  wavelength  as  that  used  in  the  recording  procedure,  the  object  image  is 
reproduced  in  the  object  plane.  Holograms  can  be  used  to  satisfy  complex  optical 
design  requirements  which  would  be  impractical  to  build  using  conventional  lens 
optics.  In  the  design  of  a  small  instrument  with  complex  optics  needed  for 
simultaneous  measurement  of  scattering  and  attenuation,  the  holographic  system  is 
attractive  since  multiple  optical  functions  can  be  incorporated  into  a  single  hologram. 

It  is  desirable  to  build  an  instrument  that  measures  ocean  optical  parameters  in  the 
blue  green  spectral  region  and  since  the  current  diode  technology  does  not  presently 
provide  compact  blue-green  lasers,  other  sources  were  investigated.  Thus,  the 
requirement  to  use  a  coherent  source  in  a  holographic  design  was  examined  early  in 
the  study.  Optical  ray  tracing  models  were  used  to  examine  the  use  of  an  incoherent 
source,  such  as  a  LED,  in  combination  with  holographic  optics.  The  results  of  these 
studies  suggested  that  this  combination  would  be  a  highly  inefficient  means  (less  than 
0. 1  %)  to  project  a  simple  beam  needed  in  a  transmissometer  application.  Under  ideal 
conditions  holograms  can  reach  efficiencies  of  ten  percent  or  greater,  but  only  with 
coherent  sources.  Presently,  inexpensive  visible  diode  lasers  are  available  only  in  the 
red  (640  and  670nm).  The  technology  is  expected  to  produce  inexpensive  blue-green 
lasers  within  three  to  five  years. 


2.3  PROPOSED  HOLOGRAPHIC  INSTRUMENT  CONCEPT 


The  proposed  design  for  the  expendable  optical  sensor  includes  four  components 
each  which  requires  specialized  holographic  design  elements.  These  are:  small  angle 
scatterometer  (0  -  10  deg.),  large  angle  scatterometer  (10  -  170  deg.),  backscatter 
meter  (170  -  180  deg.),  and  the  beam  transmissometer.  Together  these  instrument 
components  will  allow  measurement  of  the  beam  attenuation  coefficient  and  the 
volume  scattering  function.  The  development  of  holographic  elements  for  this 
instrument  permits  a  compact  design  although  the  possibility  remains  that  optical  cross 
talk  could  greatly  diminish  performance.  The  basic  four  components  of  the  proposed 
instrument  are  shown  in  Figure  2. 

2.3.1  Small  Angle  Scattering  Meter 

The  current  research  has  focussed  on  the  development  of  a  small  angle  meter  (see 
Figure  2a)  which  measures  the  MTF.  The  MTF  can  in  turn  be  used  in  combination 
with  the  beam  attenuation  measurement  to  derive  an  estimate  of  the  volume  scattering 
function  for  scattering  angles  less  than  ten  degrees.  In  this  measurement  concept  a 
bar  pattern  embedded  in  the  hologram  is  projected  through  the  sample  medium  either 
as  collimated  or  focussed  light  onto  a  linear  detector  array  which  is  used  to  measure 
the  MTF.  In  this  approach  it  is  not  necessary  to  measure  the  MTF  over  a  wide  range 
of  spatial  frequencies.  Measurements  made  at  critical  spatial  frequencies  are  used  to 
parameterize  a  semi-analytical  MTF  characteristic  of  the  sample  water  type. 

Numerical  techniques  are  used  to  invert  the  MTF  to  estimate  small  angle  scattering 
components  of  the  volume  scattering  function. 

In  the  present  instrument  concept  there  are  several  arguments  for  using  the 
focussed  versus  collimated  light.  In  the  focussed  case  the  interface  refraction  effects 
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Figure  2.  Expendable  Scatterometer  and  Transmissometer  Concept 
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are  also  recorded  on  the  same  hologram  and  therefore  are  completely  accounted  for  in 
the  projected  beam.  Construction  of  a  collimated  holographic  system  is  considered 
more  difficult  than  using  focussing  techniques.  A  hologram  which  projects  an  image 
of  a  slit  or  bar  pattern  in  prefect  collimation  along  the  sample  path  is  very  difficult  to 
construct.  Furthermore,  a  second  hologram  or  equivalent  cylindrical  lens  would  be 
required  to  focus  the  collimated  light  effectively  onto  the  linear  detector  array.  In  the 
case  of  a  focussed  light  and  holographic  optics,  all  of  the  light  is  efficiently  directed 
to  the  detector  array  as  is  the  light  scattered  in  small  angles.  Problems  of  stray  light 
are  minimized.  Other  numerical  arguments  are  made  in  the  following  sections. 

2.3.2  Large  Angle  Scattering  Meter 

It  is  envisioned  that  the  large  angle  scatterometer  (see  Figure  2b)  will  consist  of 
reflection  holograms  designed  as  long  narrow  strips  each  capable  of  focussing  the 
incident  light  received  over  a  narrow  range  of  incident  angles  to  a  single  detector. 

The  strip  holograms  can  be  laid  along  side  one  another  to  form  a  collector  bar  having 
as  much  as  three  quarters  of  the  length  of  the  sample  volume.  Multiple  bar 
reflectance  holograms  can  be  developed  on  the  same  holographic  plate.  It  is 
envisioned  that  the  bar  holograms  could  be  incorporated  into  the  open  connecting 
framework  that  surrounds  the  sample  volume. 

2.3.3  Beam  Transmissometer 

The  beam  attenuation  coefficient  is  linear  with  mass  or  volume  if  the  particle  size 
distribution  and  absorption  characteristics  are  constant.  In  a  beam  attenuation  meter  a 
point  source  of  light  is  placed  at  the  focal  point  of  a  thin  lens  forming  a  collimated 
beam  of  light.  The  beam  passes  through  the  sample  to  a  second  lens  which  focusses 


the  collimated  beam  onto  a  small  detector  placed  at  the  focal  point.  Any  light  that  is 
absorb  by  the  medium  is  lost.  Light  that  is  scattered  from  the  beam  does  not  focus  on 
the  detector.  Thus,  over  the  sample  path,  light  received  by  the  detector  measure  the 
total  beam  attenuation. 

The  design  path  length  for  making  beam  transmission  measurements  in  the  blue- 
green  region  should  be  at  least  one  meter  in  length  to  obtain  a  one  percent  accuracy 
for  the  beam  attenuation  coefficient.  Thus,  in  an  expendable  transmissometer  design 
the  collimated  beam  must  be  folded  and  passed  a  number  of  times  through  the  sample 
medium.  Our  proposed  design  is  based  on  passing  such  a  narrow  collimated  beam 
back  and  forth  between  two  holograms  each  containing  the  properly  designed  lenslet 
array.  Although  the  diagram  in  Figure  2c  shows  three  passes,  the  transmissometer 
design  will  have  a  minimum  of  four  passes  and  a  maximum  of  eight.  For  a  large 
number  of  passes  interface  losses  will  prevent  accurate  measurements.  Development 
of  a  lenslet  array  is  ideal  since  one  can  precisely  and  easily  locate  the  optics  within 
the  integrated  HOE.  In  the  design  the  primary  hologram  will  be  used  to  form  a 
collimated  beam  which  is  passed  through  the  sample  medium  to  a  second  HOE.  The 
beam  is  then  focussed  onto  a  small  plane  mirror  and  pin-hole.  The  beam  is  reflected 
to  another  lens  of  the  array  which  re-collimates  the  remaining  light  for  transmission 
along  the  return  path.  With  proper  design  of  these  optics  much  of  the  error  in  the 
attenuation  coefficient  due  to  small  angle  scattering  is  eliminated.  We  propose  to 
generate  this  holographic  pair  of  lenslet  arrays  and  demonstrate  transmissometer 
performance  in  the  laboratory. 

2.3.4  Backscatter  Meter 

The  backscatter  meter  design  is  shown  as  Figure  2d.  In  the  instrument  concept 
design  this  device  will  be  positioned  behind  the  initial  beam  forming  lens  in  the 
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transmissometer  HOE.  As  shown  in  Figure  2d  a  portion  of  the  incident  source  beam 
on  the  hologram  is  transmitted,  a  portion  reflected,  and  a  portion  transmitted  as  a 
collimated  beam.  The  light  backscattered  at  angles  near  180  degrees  will  be 
transmitted  through  the  hologram  with  approximately  50%  being  directed  back  toward 
the  source  and  50%  unaltered  in  direction.  This  latter  portion  can  be  diffracted  in  a 
second  hologram  to  separate  detectors  depending  on  the  original  scattering  angle  in 
sample  volume. 

2.3.5  Edge  Illuminated  HOEs 

The  HOEs  needed  in  backscatter  meter  and  beam  transmissometer  could  be  highly 
efficient  and  at  the  same  time  compact  by  using  edge  illuminated  HOEs.  The 
diffraction  efficiency  of  these  elements  can  be  made  nearly  100%  for  polarized  light. 
In  this  design,  light  emerging  from  the  HOE  becomes  the  collimated  beam  needed  for 
the  transmissometer.  As  above,  50%  of  the  unpolarized  backscattered  light  is 
directed  to  a  second  HOE  which  is  similar  to  the  first  but  at  right  angles  to  the  first 
so  as  to  efficiently  diffract  the  light  to  a  detector.  The  advantage  of  using  edge- 
illuminated  HOEs  is  that  they  are  monolithic  and  therefore  avoid  air-glass  interface 
problems  (Upatnieks,  1988). 
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3.0  SMALL  ANGLE  SCATTERING  THEORY 


The  fundamental  description  of  light  scattering  by  particles  in  sea  water  is  the 
volume  scattering  function,  a{d)d<j3,  which  is  defined  as  the  amount  of  light  scattered 
per  meter  into  a  differential  solid  angle  in  a  specific  direction.  For  single  scattered 
light  an  instrument  can  be  designed  to  collect  the  light  at  a  specific  angle  such  as  the 
Austin  and  Petzold  [1975]  ALSCAT  instrument.  The  problem  is  getting  reliable 
measurements  at  very  small  angles  needed  to  define  the  diffraction  peak  where  the 
physical  limits  prevent  fixed  angle  measurements.  Multiple  scattering  at  these  small 
angle  may  also  contribute  to  significant  errors.  While  measurements  made  with  the 
ALSCAT  instrument  have  been  used  extensively,  and  in  this  study,  to  estimate  the 
entire  volume  scattering  function  there  are  very  few  measurements  available  for  the 
world  oceans.  At  small  angles  the  scattering  function  has  typically  a  very  large 
forward  peak  resulting  from  (1)  diffraction  of  light  around  particles  with  diameters 
larger  than  the  wavelength  and  (2)  refraction  in  particles  with  a  nearly  water  matching 
refractive  index.  In  this  study  the  inherent  optical  parameters  refer  exclusive 
scattering  by  particles.  Thus, 


(3) 


and 


(4) 


where  b„  is  the  scattering  coefficient  for  pure  sea  water. 
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3.1  MODULATION  TRANSFER  FUNCTIONS 


In  the  case  of  significant  multiple  scattering  one  must  either  derive  the  point 
spread  function  (PSF),  (f(0,r)),  or  the  MTF  (F(^,r))  to  gain  information  on  small 
angle  scattering.  The  MTF  describes  the  loss  of  contrast  as  a  function  of  spatial 
frequency  or  the  loss  of  contrast  as  a  sinusoidal  bar  chart.  The  small  angle  scattering 
essentially  fills  in  the  vacant  spaces  between  the  individual  bars.  The  number  of 
cycles  per  unit  length  in  the  chart,  i;  relates  to  the  angular  spatial  frequency,  ^ 
(cyc/radian),  by  the  distance  r  (i.e.,  i/'=ur). 

The  PSFs  and  the  MTFs  are  integrated  by  the  two  dimensional  fourier  transform. 
Both  of  these  functions  are  circular  symmetric  and  can  be  written,  as  above,  in  terms 
of  the  distance  r  along  the  beam  and  off  axis  angle  B.  The  two  dimensional  Fourier 
transform  of  a  circular  symmetric  function  reduces  to  a  one  dimensional  Fourier- 
Bessel  integral. 

m 

f  (0,r)  =2icJj’g(27t0Tlf)F(i|»,r)ilrdt|f  (5) 

0 

F(\lr,  r)  =2nji7^,(27t0i|f)  f  (0,  r)  0d0  (6) 

0 

The  small  angle  approximation  has  been  used  in  the  second  integral. 

3.2  SMALL  ANGLE  SCATTERING  TRANSFORMS 

Since  the  PSF  is  a  function  of  the  propagation  distance  r  and  the  scattering  angle 
it  cannot  be  directly  related  to  the  volume  scattering  function.  However,  Wells 
[1969,  1972]  has  shown  a  transform  that  allows  one  to  derive  a(&)  from  the  MTF  or 
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vice  versa.  Del  Grosso  [1975]  used  this  technique  to  measure  the  MTF  for 
polydispersed  particle  size  distribution.  The  MTF  may  be  described  in  terms  of  a 
decay  function  D(^)  which  is  independent  of  measurement  distance. 


F(i|r,r)  (7) 

D(ilr)  =-r"^LoggF(t|r,r)  (8) 

The  PSF  can  result  from  two  types  of  propagation,  focussed  or  collimated  light.  It  is 
therefore  appropriate  to  think  in  terms  of  two  different  decay  functions;  one  for 
focussed  light  and  one  for  collimated  light.  The  collimated  decay  function  can 
be  directly  related  to  the  Fourier  Bessel  transform  of  a{6),  H(^). 

OD 

H(f)  =27ij<7o^2ic0!|f)  o  (0)  6d0  (9) 

0 

The  function  H(^)  is  the  MTF  of  the  scattering  function.  Although  only  values  of 
a(6)  for  ^  <  1(P  effectively  contribute  to  H(\^)  the  integration  must  be  carried  further 
to  avoid  the  problems  of  truncation  and  produce  a  satisfactory  numerical  integration. 
If  a{0)  has  been  normalized  such  that 

o  (0)  =jbpO„(0)  (10) 

and 
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l=2njo„(0)  sinSdO  (11) 

0 


then  describes  the  contributing  effects  of  scattering  to  contrast  decay  in  a  bar 

pattern  with  spatial  frequency.  At  very  large  values  of  yp  all  bar  patterns  are 
completely  obscured  by  small  angle  scattering  within  a  very  short  propagation 
distance  and  at  the  upper  limit,  yp  =  ao  and  1-H(\^')  =  1.  At  the  low  frequency  limit 
(yp  =  0)  scattering  produces  no  loss  (i.e.,  1-H„(0)=0).  The  function 


(12) 


describes  decay  of  spatial  frequencies  for  collimated  light  for  a  distribution  of 
particles  with  a  scattering  coefficient  of  bp  and  normalized  scattering  function  a^(0). 
De(^)  can  be  used  to  derive  the  normalized  volume  scattering  function. 


a„(0)  =2itJt7o(27i0;|f)H„(i|i)i|rdiji 

O 


(13) 


For  the  case  of  a  focussed  spatial  decay  function  the  Dc(^)  function  is 
integrated  over  the  [0,^]  interval  with  the  dummy  variable  u  to  give  a  larger  value  of 
decay  demonstrating  the  reduced  effect  of  scatters  near  the  focal  point 
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T 

D^(i|r)  du 


The  reverse  transform  is  simply 


_  (i|f)  ] 


One  can  also  think  of  an  equivalent  relationship  in  the  6  domain.  Since  the 
Fourier-Bessel  transform  of  d[i^H(i/')]/d0  is  -d[(?<r(0)]/d0  one  can  derive  the  effective 
scattering  function  for  focussed  light,  From  the  above  relationships  it  is 
apparent  that 


_  d[ilr/ff  (i|r)  ] 


d[0af(0)]  “r 

- ^ - =-2Tt|i7„{27i0l|r)Hp(»lr)i|fdi|j 


implying 


d0 
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and 


0^(0)  =0'^Jo^(u)  du  (19) 

e 

The  physical  interpretation  of  this  result  is  simply  for  the  collimated  case  light 
that  is  scattered  at  the  same  angle  along  the  beam  will  lie  at  the  same  position  in  the 
PSF  fc(0,r)  whereas  in  the  focussed  case  light  scattered  at  the  same  angle  will  end  up 
in  different  locations  or  conversely  light  received  at  a  single  position  on  will 
have  come  from  angle  6  and  greater  angles.  Thus,  a0)  represents  the  average  value 
of  a^B)  over  the  interval  (0,ir). 

Thus  an  alternative  and  equivalent  transform  for  the  focussed  case  gives 
from  OfiB)  which  is  obtained  form  0^(B)  using  equation  (18).  The  focussed  decay 
transform  for  becomes 

0* 

=Cp-2it|jr^(27t0i|f)Of(0)0de  (20) 


and  the  inverse  transform 


a^(d)  =2nJj^(2ndi^)  [£>^(i|r) -Cp]  i|fdi|f 

O 


(21) 


Use  of  the  latter  equation  will  tend  to  reduce  noise  in  estimating  the  0f(B)  since 
one  avoids  introduction  of  numerical  noise  in  the  differentiation  step  needed  in  the 
earlier  transform. 
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3.3  COMPUTATIONAL  TECHNIQUES 


The  small  angle  scattering  analysis  procedures  are  shown  in  the  block  diagram  of 
Figure  3.  Three  basic  analysis  procedures  were  followed  leading  to  a  derived  small 
angle  scattering  function.  (1)  In  the  laboratory  experiments  data  were  collected  on  a 
linear  detector  array  and  used  to  estimate  contrast  loss  at  selected  spatial  frequencies. 
These  measurement  data  were  then  used  to  estimate  the  medium  MTF  and  finally 
information  on  small  angle  scattering.  (2)  A  parallel  analysis  was  conducted  starting 
with  the  particle  size  for  the  laboratory  experiment  and  using  Mie  theory  to  estimate 
the  volume  scattering  function  (VSF).  The  VSF  was  used  in  an  instrument  model  to 
simulate  the  laboratory  measurements.  The  VSF  was  transformed  to  the  MTF.  The 
estimated  MTF  was  then  compared  to  the  laboratory  measurements  and  simulated 
measurements.  (3)  A  VSF  based  upon  measurement  data  was  used  to  calculate  the 
MTF.  Instrument  size  parameters  were  then  used  in  the  model  to  simulate  measured 
contrast.  These  measurements  were  then  compared  to  the  MTF  derived  from  VSF 
and  furthermore  were  used  to  derive  an  estimate  of  the  MTF.  The  inverse  transform 
was  then  applied  to  derive  the  small  angle  portion  of  the  VSF.  Three  basic  computer 
codes  were  developed  to  simulate  the  processes  for  a  holographic  small  angle 
scattering  meter.  The  combination  of  codes  allows  end  to  end  measurement  analysis 
of  the  volume  scattering  function. 

A  code  was  developed  to  make  the  Fourier-Bessel  integration  transform 
calculations  described  above.  This  code,  MTF_SCAT,  uses  an  8000  point  Labatto 
quadrature  in  the  Fourier-Bessel  subroutine  and  cubic  splines  in  estimating  the  local 
derivative.  This  procedure  provides  both  forward  (VSF  to  MTF)  and  inverse 
calculations  (MTF  to  VSF).  Equi-space  sampling  in  the  log  domain  was  selected 
since  both  functions  are  strongly  peak  at  the  origin.  The  VSF  developed  for  clear 
ocean  water  and  used  in  these  analyses  was  based  upon  Petzold’s  [1972]  scattering 
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Figure  3.  Small  Angle  Scattering  Analysis  Procedures 
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data  at  530nm  collected  for  the  Tongue  of  the  Ocean.  The  measured  value  of 
attenuation  for  this  sanpiple  was  0. 199m  ‘  and  the  total  scattering  coefficient  was 
estimated  to  be  0.177m  ‘.  These  data  were  fitted  [Tanis  et.  al.,  1986]  to  a  smooth 
VSF,  as  shown  in  Figure  4,  using  Pede  splines.  In  the  MTF  calculations  the  VSF 
was  retained  for  angles  less  than  100°.  At  larger  angles  the  VSF  was  tapered  to  zero 
and  the  entire  modified  function  was  renormalized.  This  procedure  provides  a  smooth 
function  necessary  for  the  Fourier-Bessel  transform  calculations.  While  it  is  realized 
that  only  forward  scattering  (d<  10°)  is  important  to  the  observed  MTF  a  smooth 
truncation  must  be  provided  to  minimize  truncation  effects  which  are  more  apparent 
in  the  inverse  calculations.  The  derived  MTF  for  the  AUTEC  sample  is  shown  in 
Figure  5. 

The  second  code  developed,  PROBE,  is  a  Monte  Carlo  model  with  fixed  spatial 
frequencies  designed  to  simulate  the  measurement  instrument  by  projecting  a  bar 
pattern  through  the  scattering  medium  and  focussing  the  pattern  modified  by  small 
angle  scattering  on  to  a  linear  detector  array.  The  approach  taken  was  to  force 
collisions  within  the  sample  medium  to  avoid  counting  those  which  escaped.  The 
fraction  of  photons  escaping  the  first  collision  and  passing  through  the  medium  is  e"^, 
where  t,  is  the  optical  thickness  of  the  medium.  One  forces  a  collision  by  sampling 
from  the  density  function. 


p(x)  dx  = 


e''dx 


(22) 


The  potential  bias  is  handled  by  multiplying  a  weight  by  1-e'",  which  yields  the 
fraction  of  photons  remaining  in  the  packet.  The  weight  is  reduced  at  each  collision 
until  it  falls  below  a  cutoff.  At  this  time  a  new  photon  is  started  in  the  simulation. 
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Contrast 


Theta  (radians) 


- Modified  VSF - Orig  VSF 


Figure  4.  Modified  Volume  Scattering  Function  for  AUTEC  Sample 
Number  7  and  Original  Fitted  Curve  Using  Petzold's 
Measurement  Data  [1972] 


Derived  MTFI  - alpha=  .34421 


Figure  5.  Modulation  Transfer  Function  for  AUTEC  Sample  With  Path 
Length  of  25  cm.  Also  Shown  is  the  Fitted  MTF  (Solid  Curve) 

With  a  =  0.3442,  p  =  .2428,  ;  =  0.39,  and  y  =  1 .00 
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Model  scattering  includes  both  Rayleigh  and  Mie  although  for  the  present  study  only 
Mie  calculations  were  made. 

The  spatial  distribution  of  photons  leaving  the  hologram  was  uniform  since  every 
position  in  the  hologram  contributes  light  to  the  image.  The  selection  of  photon 
initial  direction  was  based  upon  the  distribution  of  unscattered  photons  in  the  image 
plane  as  calculated  from  the  sinusoidal  bar  pattern  (see  Figure  6).  The  bar  pattern 
selected  was  consistent  with  that  chosen  for  the  laboratory  experiments  and  consists  of 
1024  detector  elements  (as  in  the  actual  linear  array).  The  extensions  on  either  side 
of  the  central  pattern  provide  calibration  in  the  contrast  calculations  and  help  to 
prevent  humping  of  the  scattered  portion  due  to  apparent  unequal  lateral  scattering. 

For  a  photon  focussed  on  the  end  detector  only  energy  scattered  in  one  direction 
along  the  array  will  be  recorded  whereas  in  the  center  recording  will  be  made  in  both 
directions. 

In  the  25cm  instrument  simulations  the  hologram  was  4.0  x  4.0cm,  the  detector 
array  size  at  image  focus  was  0.04  X  10.24cm.  It  was  assumed  that  optics  were 
present,  as  in  the  laboratory  experiments,  to  re-focus  the  bar  pattern  onto  an  actual 
array  (.01  X  2.56cm).  These  latter  optics  were  not  part  of  the  present  simulations. 
Examples  of  the  scattered  portion  for  a  25cm  path  length  instrument  using  the 
AUTEC  No. 7  volume  scattering  function  are  shown  in  Figure  7  for  sample  runs  of 
200,  2500,  and  7500  photons.  Near  convergence  is  achieved  at  Np  =2500  for 
AUTEC  No.7  and  is  sufficient  to  study  the  behavior  of  scattered  energy  in  the  present 
study.  More  rapid  convergence  was  observed  for  other  volume  scattering  functions 
with  a  smaller  diffraction  peak.  The  simulator  is  completely  parameterize  to  allow 
the  user  to  set  up  any  instrument  geometry,  projected  pattern,  or  scattering 
conditions. 

A  third  code  contained  three  procedures;  (1)  to  derive  a  semi-analytic  fit  to  the 
calculated  MTF,  (2)  to  generate  a  series  of  MTF’s  based  upon  a  single  parameter 
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Figure  6.  Simulated  Unscattered  Energy  on  the  Linear  Detector  Array  for 
a  Path  Length  of  25  cm  and  Beam  Attenuation  of  0.2  m-i 


Figure  7.  Comparison  of  Simulated  Scattered  Energy  on  Linear  Detector 
Array  Using  200,  2500,  and  7500  Photons 
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input,  a,  and  (3)  to  calculate  the  best  fit  based  upon  measurement  data.  The  MTF 
curve  fitting  technique  involves  fitting  two  exponential  functions  which  are  then 
splined  together.  The  first  function  fi(^)  is  a  fit  to  the  MTF  at  low  spatial 
frequencies. 


(i|r)  =  (1-e*^) 


(23) 


where  t  is  the  product  of  path  length  and  c,  rp  is  spatial  frequency  in  cycles  per 
radian,  and  7  is  a  smoothing  constant.  The  function  f,  equals  1  at  (^^=0)  and  e"^  at 
large  spatial  frequencies.  The  second  function  has  the  same  formulation  as  the 
first  and  describes  the  MTF  at  large  spatial  frequencies. 

f2(i|r)  (24) 

where  c  is  a  smoothing  constant  for  the  second  function.  The  function  f2  converges  to 
for  large  values  and  1  equals  at  =  0. 

Initially  both  functions  are  estimated  using  a  known  t  and  a  MTF  characteristic  of 
the  water  type,  i.e,  one  must  assume  a  water  type  and  scattering  function  to  get  a 
characteristic  MTF.  The  two  functions  are  combined  into  a  smooth  composite 
function  with  spline  functions.  An  example  fit  is  shown  in  Figure  5  for  an  a  = 
0.3442  and  /3  =  .2428.  An  RMS  error  of  5.1X10“*  was  calculated  for  this  fit.  Either 
function  can  be  changed  slightly  and  the  two  functions  then  reconnected.  For  MTF 
fitting  routine  the  fi  function  is  assumed  constant  for  spatial  frequencies  less  than  a 
few  tenths  of  a  cycle/radian.  The  f2  function  is  allowed  to  vary  using  only  the  a 
parameter.  Example  MTFs  are  shown  in  Figure  8  for  the  AUTEC  sample  and  a  t 


®>ERIM  _ 

value  of  1.  These  three  MTF’s  were  then  inverted  as  shown  in  Figure  9  to  indicate 
the  corresponding  variability  in  small  angle  VSF  at  angles  less  than  0.4  degrees. 

Also  shown  is  the  original  AUTEC  VSF  which  compares  well  with  the  inverted 
functions  at  large  angles.  The  noise  observed  at  angles  greater  than  one  radian  is  an 
artifact  of  the  Bessel  function  transform. 
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Spatial  Frequency  (cyc/rad) 


-  alpha=  0.4511  - alpha=  0.4060 - alpha=  0.4962 


Figure  8.  Modulation  Transfer  Functions  Obtained  for  Selected  Values 
of  Alpha 


- AulecNo.7  - alpha=  0.4511  — —  alpria=  0.4060  - alpha=  0,4962 


Figure  9.  Volume  Scattering  Functions  Derived  for  the  Three  MTF  s 
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4.0  LABORATORY  METHODS 


An  experimental  setup  was  assembled  to  use  both  methods  in  the  recording  of  the 
hologram  and  also  for  pattern  projection.  The  laboratory  laser  table  setups  are  shown 
as  Figures  10  and  1 1  respectively.  The  primary  laser  was  a  Spectra  Physics  125 
HeNe  laser  with  90  to  lOOmW  output  and  was  used  to  record  and  play  back  the 
hologram  (see  Figure  12).  A  second  lOmW  laser  was  used  to  measure  particle 
attenuation  in  the  sample  tank  (see  Figure  13). 

4.1  HOLOGRAM  RECORDING  SETUP 

The  hologram  recording  material  was  2”X2'’  Kodak  Spectroscopic  plates  with 
photo  sensitive  backing.  As  shown  in  Figure  10  the  beam  was  split  with  the  aid  of  a 
variable  spitter  into  a  beam  to  illuminate  the  bar  transparency  through  the  water  filled 
tank  and  a  second  to  create  the  interference  pattern  on  the  recording  plate.  The 
plexiglass  tank  had  inside  dimensions  of  25  by  100  cm  and  was  constructed  of  3/8" 
clear  plexiglass.  Scattering  within  the  plexiglass  sheet  was  not  considered  a  problem. 
For  hologram  recording  the  tank  was  filled  with  nearly  50  liters  of  filtered 
deionized/distilled  water  made  available  from  ERIM’s  Electronics  Facility. 

4.2  LABORATORY  OPTICAL  SYSTEM  FOR  MEASUREMENTS 

In  the  process  of  recording  the  sinusoidal  bar  pattern  all  surface  refraction  effects 
are  included.  A  field  lens  was  used  to  gather  the  light  at  the  image  plane  and  pass  it 
to  the  linear  array  via  an  imaging  lens.  The  linear  detector  array  was  an  EG&G 
model  RL  1024K  unit  with  1024  element  with  center  spacing  of  25  microns  and  a 
length  of  100  microns.  The  bar  pattern  and  optics  were  sized  to  efficiently  image  the 
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Figure  10.  Optical  System  for  Hologram  Recording 


Figure  1 1 .  Optical  System  for  Light  Scattering  Measurement 
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pattern  through  the  tank  and  onto  the  linear  array.  The  detected  energy  on  the  linear 
array  was  sampled  using  an  array  dump  routine  placed  on  the  DSP56  signal 
processing  board.  The  array  data  were  digitized  and  passed  to  a  286  PC  for  display 
and  data  analysis.  A  flow  chart  of  the  data  recording  system  is  shown  in  Figure  14. 
PC  software  (LINEAR IK)  was  developed  to  provided  display  of  the  scan  profile  at 
two  different  magnifications.  Data  were  recorded  with  sixteen  bits  of  dynamic  range. 
Electronic  noise  alone  reduced  measurement  significance  to  about  fourteen  data  bits. 
Each  array  scan  integrated  the  detected  light  energy  over  a  10  ms  period.  Multiple 
scans  were  then  averaged  in  the  PC  data  acquisition  software  to  reduce  electronic  and 
other  systematic  sources  of  noise.  When  data  profiles  were  averaged  over  400  scans 
about  1 1  data  bits  were  found  to  be  significant.  This  precision  was  considered  to  be 
adequate  to  record  the  dynamic  range  of  the  projected  bar  pattern.  Some  drift  in  the 
system  was  observed  so  that  noise  reduction  by  averaging  had  little  consequence  for 
collections  involving  more  than  400  scans. 

4.3  SINUSOIDAL  PATTERN 

The  sinusoidal  bar  patterns  were  computer  generated  and  transferred  to  film  using 
a  high  resolution  image  filming  recorder.  The  spatial  frequencies  selected  were 
designed  for  the  25  cm  path  length  and  to  characterize  the  portion  of  the  MTF 
sensitive  to  small  angle  scattering.  The  spatial  frequencies  of  this  pattern  depend  on 
the  distance  between  hologram  and  image  plane.  The  path  between  the  hologram  and 
image  plane,  in  the  laboratory,  was  found  to  be  41  cm  including  the  25  cm  water  path 
and  two  plexiglas  sheets.  The  resulting  spatial  frequencies  are  given  in  Table  1  with 
units  of  cycles  per  radian. 
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Figure  14,  Laboratory  Data  Recording  System 
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Table  1. 

Bar  Pattern  Spatial  Frequen^'ies 
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Frequency  (cyc/rad) 
25.6 
51.3 

102.5 
205.0 

341.5 


Number  of  Cycles 
1 
2 
2 
3 
5 


4.4  SCATTERING  MEDIUM 


The  tank  was  seeded  with  polystyrene  spheres  available  from  Polyscience  Inc. 
under  the  trade  name  POLYBEAD.  The  concentration  of  spheres  was  diluted  into 
two  fractions  with  concentrations  of  1.1  X  10’  and  1.1  X  10*  particles  per  milliter. 
Aliquots  from  these  dilutions  were  then  transferred  to  the  tank  which  contained  47.5 
liters  of  water.  The  spheres  are  considered  non-absorbing  and  have  a  nominal 
refractive  index  of  1.60-fOi. 

4.5  HOLOGRAM  SNR  AND  DYNAMIC  RANGE 


Several  holograms  were  developed  for  this  experiment  with  the  objective  being  to 
create  a  hologram  with  high  efficiency  and  low  signal  to  noise  ratio  (SNR). 
Performance  characteristics  of  several  of  these  holograms  were  measured  in  the 
laboratory.  Minimum  SNR  was  p)erhaps  of  greater  impc .lance  than  maximum 
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efficiency  since  an  early  decision  was  made  to  use  the  large  HeNe  laser  for  hologram 
construction. 

The  dynamic  range  of  the  input  bar  pattern  recorded  on  the  hologram  was  10*:  1. 

It  was  expected  that  the  resulting  hologram  would  produce  a  dynamic  signal  range,  in 
the  projected  image,  in  excess  of  10^:1.  This  range  was  needed  to  accurately  measure 
changes  in  contrast  between  the  light  and  dark  fields  due  to  particle  scattering.  The 
best  measured  dynamic  range  between  the  light  and  dark  image  fields  without  the  tank 
and  scattering  medium  was  600: 1 .  When  the  hologram  was  mounted  for  playback 
through  the  tank  this  number  was  reduced  to  100:1.  The  former  range  is  indicative  of 
the  inherent  SNR  of  the  hologram  due  to  speckle.  The  speckle  size  was  estimated  to 
be  6.6  microns  square  at  the  HeNe  633nm  wavelength.  The  number  of  speckles 
located  in  one  detector  element  at  the  image  plane  was  estimated  to  be  2.3X10^. 

Thus,  a  measure  of  the  speckle  SNR  of  the  hologram  is  approximately  48:1  which  is 
indicative  of  the  signal  standard  deviation  for  a  uniform  light  field.  The  600: 1  ratio  is 
indicative  of  the  dynamic  range  of  the  hologram  where  the  minimum  is  due  largely  to 
speckle.  The  reduction  in  this  ratio  through  the  tank  and  water  medium  is  due  to  both 
scattering  in  the  water  and  reception  of  stray  light  into  the  detector  field.  The  speckle 
noise  is  a  feature  of  the  hologram  itself  so  that  if  the  structure  remains  fixed  the  noise 
limit  for  a  large  number  of  samples  will  approach  some  other  smaller  limit.  For  the 
present  experiment  this  lower  limit  is  determined  by  the  number  of  significant  bits  in 
the  resulting  measurements  (approximately  11).  The  signal  due  to  scattering  and  its 
modulation  by  the  sinusoidal  pattern  were  estimated  to  be  1%  to  10%  of  the  dark 
field  signal.  With  careful  control  of  stray  light  and  the  speckle  problem,  SNRs  of 
10^:1  should  be  achievable. 
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5.0  EXPERIMENTAL  AND  SIMULATION  RESULTS 

The  laboratory  experiments  were  intended  to  demonstrate  the  holographic 
technique  and  cannot  be  considered  quantitative.  The  instrument  scattering  model  was 
used  for  two  purposes.  First,  to  simulate  contrast  loss  in  a  sinusoidal  pattern  for 
specific  monodispersions  of  polystyrene  spheres  used  in  the  laboratory  experiments 
and  second,  to  demonstrate  capability  of  the  proposed  25cm  instrument  to  measure 
small  angle  scattering.  In  the  second  effort  contrast  loss  was  simulated  from  the 
model  knowing  the  scattering  phase  function  for  a  polydispersion  of  ocean  particles. 
These  contrasts  were  then  used  to  estimate  a  contrast  decay  function  at  selected  spatial 
frequencies.  A  numerical  minimum  square  error  procedure  was  used  to  estimate  the 
MTF  based  upon  these  selected  spatial  frequencies.  The  derived  MTF  was  then 
inverted  using  Fourier-Bessel  transforms  to  estimate  the  small  angle  portion  of  the 
volume  scattering  function. 

5.1  POLYSTYRENE  SPHERE  SEEDING  EXPERIMENTS 

A  seeding  experiments  were  conducted  for  four  particle  diameters:  0.049^m, 
0.99/im,  2.92/im  and  25.0/im.  In  each  of  these  experiments  polystyrene  spheres  were 
mixed  with  ERIM  distilled  water  in  five  or  six  steps  resulting  in  a  wide  range  of 
panicle  concentrations.  A  average  scan  profile  recorded  on  the  detector  array  is 
shown  in  Figure  15.  All  but  the  first  cycle  of  the  pattern  is  shown  which  could  not 
be  included  in  the  focussed  image.  The  reduced  peak  size  at  higher  frequencies  is  an 
artifact  of  the  photographic  and  hologram  recording  processes  and  could  not  be 
avoided.  The  dynamic  range  of  the  detector  was  maximized  for  these  observations 
and  was  far  in  excess  of  hologram  and  stray  light  noise. 
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The  average  background  attenuation  at  633nm  was  estimated  to  be  0.320  m-1. 
The  accepted  value  at  665nm  is  0.364  m-1  [Zaneveld,1984].  The  range  in  the 
attenuation  coefficient  for  these  experiments  ranged  from  0.32  to  1.0  m  *.  Beam 
attenuation  measurements  for  the  2.92  micron  polystyrene  spheres  is  shown  in 
Figure  16. 

Calculated  contrasts  are  plotted  in  Figures  17  and  18  for  the  observed  spatial 
frequencies  and  attenuation  values.  The  25/im  diameter  spheres  showed  the  largest 
reduction  in  contrast  of  approximately  0.02  from  the  smallest  to  largest  frequency. 
For  the  0.49/im  spheres  contrast  reduction  was  minimal.  The  results  from  the 
0.99/im  diameter  experiments  were  noisy  and  inconclusive.  The  measurement  series 
with  the  2.92/xm  spheres  produced  small  contrast  reductions  on  the  order  of  0.005. 
These  results  are  not  directly  comparable  to  the  simulated  results  but  do  show  that 
some  contrast  reduction  was  observed  with  the  larger  sphere  diameters. 


5.2  SIMULATION  OF  SEEDING  EXPERIMENTS 


A  series  of  simulations  were  made  to  aid  interpretation  of  the  laboratory 
measurements.  A  volume  scattering  function  was  calculated  for  each  particle 
diameter  using  a  Mie  code  developed  by  Kattawar  [1988]  and  a  index  of  refraction  of 
1.60+0i.  Calculated  VSFs  are  shown  in  Figure  19  for  six  particle  sizes.  As 
expected  the  small  particle  sizes  exhibit  a  broad  peak.  The  larger  the  particle  size  the 
broader  the  small  angle  peak.  Each  of  these  VSFs  was  used  to  obtain  a  simulated 
trace  of  scattered  photon  energy  on  the  linear  detector  array  using  the  laboratory 
geometry.  These  simulations  are  shown  in  Figure  20  and  indicate  that  at  the  fixed 
spatial  frequencies  of  the  experiment  the  contrast  has  decayed  to  the  high  frequency 
limit  [e'^  except  at  the  larger  particle  size  diameters.  Thus,  the  simulation  results  are 
consistent  with  the  laboratory  observations. 
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Detector  Number 


Figure  1 5.  Scan  of  Linear  Detector  Array  in  the  Seeding  Experiments  With 
2.92  Micron  Particles 


Figure  1 6.  Beam  Attenuation  Measurements  Versus  Number  of  2.92 
Micron  Polystyrene  Particles 
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jSSSa  C=0.32m  1  EB  c=0.35ni  1  ^  C=0.39m-1 
C=0.46m-1  C=0.61m-1 


(a) 


0.75 


1.3  1.3  2.5  2.5  5  5  5  8.3  8.3  8.3  8.3  8.3 

Spatial  Frequency  (cycles/cm) 


Bssa  C=0.32m-1  IM  C=0.37in-1  C=0.42m-1 

FSffl  C=0.45m-1  C=0.53in-1 

(b) 

Figure  18.  Estimated  Contrasts  (Difference/Sum)  at  Spatial  Frequencies 
Recorded  for  the  (a)  2.92  Micron  Spheres  and  (  b)  25.0  Micron 
Spheres 
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Polystyrene  Spheres  D=0.49um  Polystyrene  Spheres  D=0.99um 


Linear  Array  Detector  Number  Linear  Array  Detector  Number 


Polystyrene  Spheres  D=2.92um  Polystyrene  Spheres  D=6.44um 


Linear  Array  Detector  Number 


Unear  Array  Detector  Number 


Polystyrene  Spheres  D= 1  S.Oum  Polystyrene  Spheres  D=25.0um 


Linear  Array  Detector  Number  Linear  Array  Dacaclor  Number 


Figure  20.  Simulated  Scattered  Energy  Received  on  Linear  Detector  Array 
for  a  Path  Length  of  41  cm  and  a  Beam  Attenuation  of  0.2  m-i 
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5.3  SIMULATION  OF  AUTEC  WATER 

The  performance  of  the  holographic  small  angle  scattering  meter  was  investigated 
for  natural  waters  using  the  sensor  simulator  and  the  AUTEC  volume  scattering 
function.  The  25cm  path  instrument  was  simulated  for  each  of  two  scattering 
functions  shown  in  Figure  21  where  the  VSFl  is  the  original  scattering  function  and 
VSF2  is  a  modified  function  with  small  angle  scattering  reduced  by  approximately  a 
factor  of  ten.  The  derived  MTF  for  each  of  these  functions  is  shown  in  Figure  22 
along  with  the  simulated  measurements.  The  simulated  patterns  of  scattered  energy 
on  the  detector  array  for  VSFl  and  VSF2  from  which  contrasts  were  calculated  is 
shown  in  Figure  23a.  The  MTFl  function  was  previously  parameterized  (see  Figure 
5).  This  parameterization  was  used  to  refit  MTFl  based  upon  the  simulated 
measurement  data.  In  this  procedure  the  previously  derived  f,(^)  is  held  constant 
while  the  a  parameter  in  f2(i/')  is  allowed  to  vary  until  a  least  squares  fit  is  obtained 
with  the  measurement  data.  An  a  value  of  0.3C16  was  obtain  for  the  fitted  function 
shown  in  Figure  22.  The  simulated  measurements  made  using  VSF2  are  shown  in 
Figure  24.  The  MTFl  function  and  parameterization,  for  this  demonstration,  can  be 
considered  to  be  the  standard.  The  simulated  measurements  made  with  VSF2  were 
then  used  in  combination  with  the  standard  MTF  parameterization  to  derive  an 
estimate  of  MTF2  also  shown  in  Figure  24.  The  derived  MTF  appears  to  be  a  good 
fit  to  the  measurement  points  which  lie  slightly  above  the  MTF2  curve.  In  the  last 
step  this  derived  MTF  was  inverted  to  estimate  the  small  angle  scattering  function  for 
VSF2  as  shown  in  Figure  25.  The  resulting  scattering  function  lies  slightly  above 
VSF2  as  would  be  expected  based  upon  the  relationships  between  the  MTFs  given 
previously  in  Figure  24.  In  another  exercise  an  attempt  was  made  to  derive  the  MTF 
for  the  25.0  micron  particles  and  compare  the  results  to  simulated  and  laboratory 
measurements.  The  MTF  of  the  original  VSF  for  these  particles  could  not  be 
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-  VSF1 - VSF2 


Figure  21 .  Volume  Scattering  Function  for  AUTEC  Sample  Number  7 
(VSF1)  and  a  Modified  Version  With  Small  Angle  Scattering 
Reduced  by  Approximately  1 0x  (VSF2) 


- Derived  MTF1  - Derived  MTF2  - alpha=  .36163  —  Simulated  Data 


Figure  22.  Derived  Modulation  Transfer  Functions  for  AUTEC  Number  7 
(MTF1 )  and  for  VSF2  (MTF2)  Using  a  25  cm  Instrument 
Sample  Path  Length  and  C=0.2  m-i .  The  Dotted  Curve  With 
a  =  ,3616  was  Fitted  to  MTF1  Using  the  Simulated 
Measurements 
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(a) 


Figure  23.  Simulated  Distribution  of  Scattered  Energy  on  the  Linear 
Detector  Array  Using  Volume  Scattering  Functions  VSF1 
and  VSF2.  Panel  (a),  Path  Length  =  25  cm,  C  =  0.2  m-i  Panel 
(b)  Path  Length  =  1 00  cm,  C  =  0.2  m-i 
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Spatial  Frequency  (cyc/rad) 


-  Derived  MTFl  -  Derived  MTF2  -  alpha=  .45404  +  Simulated  Data 


Figure  24.  Derived  MTF  (Dotted  Curve)  Based  Upon  Simulated 
Measurements  Made  Using  VSF2 


- VSFt  - VSF2  - Derived  VSF2 


Figure  25.  Volume  Scattering  Functions  for  MTFl ,  MTF2,  and  Derived 
VSF  Based  Upon  the  Simulated  Measurements  Using  VSF2 
and  the  MTF1  Parameterization 
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calculated  because  of  the  large  lobes  in  the  scattering  function.  As  an  alternative  this 
VSF  was  fitted  with  a  smooth  curve  as  shown  in  Figure  26  suitable  for  MTF 
conversion.  The  resulting  MTF  is  plotted  in  Figure  27  together  with  both  the 
laboratory  measurements  and  the  simulated  measurements.  The  laboratory 
measurements  are  uncalibrated  and  therefore  are  not  directly  comparable  to  the 
simulated  data  but  do  appear  to  have  a  comparable  decay  relationship  with  spatial 
frequency. 

5.4  DISCUSSION  OF  RESULTS 

The  results  from  the  experiments  and  simulations  are  encouraging  and  support  the 
holographic  MTF  measurement  concept  for  small  angle  scattering.  Obviously  more 
quantitative  measurements  are  needed  and  more  complete  simulations  made  which 
incorporate  detector  noise  and  instrument  drift.  These  considerations  are  especially 
important  to  simulating  instrument  performance  tradeoffs  with  path  length.  In  the 
present  simulations  nearly  perfect  resolution  of  scattered  plus  direct  energy  was 
assumed.  The  simulated  pattern  of  scattered  energy  on  the  detector  array  shown  in 
Figure  23b  appear  noisy  because  the  simulation  times  were  insufficient  to  sample  all 
the  multiple  scattering  effects  with  the  longer  one  meter  path  length.  Patterns  of 
scattered  energy  for  25cm  and  lOOcm  path  lengths  otherwise  appear  quite  similar  and 
the  demonstrated  measurement  capability  above  could  be  expected  to  be  comparable 
independent  of  path  length.  For  the  present  simulations  for  Tongue  of  the  Ocean 
waters  the  ratio  of  scattered  energy  to  direct  in  a  25cm  instrument  with  a  1024 
element  array  is  about  1:500.  In  a  one  meter  instrument  this  ratio  becomes  1:95  and 
in  a  four  meter  instrument  1:16.  This  ratio  also  represents  a  reasonable  lower  limit 
for  the  SNR  resulting  from  a  combination  of  hologram  and  detector  noise.  In  the 
present  laboratory  demonstration  a  SNR  of  100:1  was  achieved  and  extensive 
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- Fitted  - Original 


Figure  26.  Volume  Scattering  Function  for  25.  0  |im  Particles  and  Fitted 
Smooth  Curve 


- Derived  MTF  □  Simulation  Data  Measured  Data 


Figure  27.  Comparison  of  Derived  MTF  From  the  Fitted  Function  for 
25.0  pm  Particles  With  Simulated  and  Laboratory 
Measurements 
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averaging  was  therefore  required  to  reduce  the  effects  of  noise.  The  hologram 
speckle  noise  is  likely  to  be  limiting  in  a  practical  instrument.  However,  a  hologram 
SNR  between  500:1  and  1000:1  is  achievable. 
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6.0  CONCLUSIONS  AND  RECOMMENDATIONS 


A  preliminary  design  concept  has  been  developed  for  a  holographic  instrument  to 
measure  scattering  and  transmission  optical  properties  of  ocean  waters.  Use  of 
holographic  optics  allows  compact  integration  needed  in  an  expendable  unit.  The 
laboratory  experiments  and  simulations  demonstrate  that  a  holographic  MTF 
instrument  is  feasible  and  will  provide  valuable  small  angle  measurements. 

It  is  recommended  that  this  basic  research  continue  and  be  directed  towards  the 
following  goals. 

1 .  Develop  methodology  to  construct  a  lenslet  hologram  needed  for  a  multiple 

path  transmissiometer  and  to  integrate  with  the  hologram  needed  for  the 

MTF  small  angle  meter. 

2.  Develop  recording  techniques  to  construct  needed  holograms  with 
SNR  values  in  excess  of  500:1. 

3.  Develop  complete  calibrations  techniques  for  the  MTF  meter  and 
conduct  laboratory  measurements  with  poly-dispersions 
characteristic  of  oceanic  waters. 

4.  Develop  simulations  for  the  transmissometer,  large  angle 
scatterometer,  and  backscatter  meter.  Compare  results  with 
theoretical  models  and  measurement  data. 

5.  Develop  a  complete  integrated  concept  design  which  includes  all 
measurement  functions. 

After  completion  of  the  concept  design  phase  the  project  can  be  transferred  to 
a  Phase  I  prototyp)e  design  which  will  involve  fabrication  of  a  retrievable  instrument. 
Laboratory  and  tank  tests  will  be  conducted  to  evaluate  accuracy  and  stability.  Phase 
II,  prototype  evaluation  will  involve  participation  in  a  series  of  Navy  ocean  optics 
cruises.  During  this  phase  an  advanced  prototype  instrument  will  be  developed  which 
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includes  the  XBT  wire  launch  technology.  During  Phase  III,  advanced  prototype 
development,  multiple  retrievable  units  will  be  develop  and  an  extensive  evaluation 
program  will  be  conducted. 
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